Astronomy & Astrophysics manuscript no. Mengel ' astroph 


February 2, 2008 


(DOI: will be inserted by hand later) 





Star-formation in NGC 4038/4039 from broad- and narrow band 
photometry: Cluster Destruction?* 

Sabine Mengel 1 , Matthew D. Lehnert 2 , Niranjan Thatte 3 , and Reinhard Genzel 2 ' 4 



in 
o 
o 



o 



> 

in 
^t- 
^j- 
»n 
o 
in 
o 

^ : 

0,; 

6 ' 
a: 



1 European Southern Observatory, Karl-Schwarzschild-Str. 2, D-85748 Garching, Germany 

2 Max-Planck-Institut fiir extraterrestrische Physik, GiessenbachstraBe, D-85748 Garching, Germany 

3 University of Oxford, Dept. of Astrophysics, Denys Wilkinson Building, Keble Road, GB-Oxford OX1 3RH 

4 Also at Department of Physics, University of California at Berkeley, 366 Le Conte Hall, Berkeley, CA 94720-7300 
email: smengel@eso.org, mlehnert@mpe.mpg.de, thatte@astro.ox.ac.uk, genzel@mpe.mpg.de 

Received ; accepted 



Abstract. An abstract should be given 



Accurately determining the star formation history in NGC 
4038/4039 - "The Antennae" is hampered by variable and 
sometimes substantial extinction. We therefore used near in- 
frared broad- and narrow-band images obtained with ISAAC 
at the VLT and with SOFI at the NTT to determine the recent 
star formation history in this prototypical merger. In combina- 
tion with archival HST data, we determined ages, extinction 
and other parameters for single star clusters, and properties of 
the cluster population as a whole. 

About 70% of the K s -band detected star clusters with 
masses >10 5 Mq are younger than 10 Myrs (this is approxi- 
mately an e-folding time for cluster ages), which we interpret 
as evidence for rapid dissolution but not free expansion. 

The total mass of K-band selected clusters is about 5 to 
lOxlO 8 Mq and represents about 3-6% of the total molecular 
gas. However, this takes into account only the detected clus- 
ters and in view of the rapid dissolution means that this is only 
a lower limit to the total mass of stars produced in clusters 
during the burst. Studies of cluster formation in other galax- 
ies recently suggested short cluster dissolution timescales, too, 
which means that star formation rates may have been severely 
underestimated in the past. 

Extinction is strongly variable and very high in some re- 
gions, but around Ay=1.3 mag on average. Even though most 
clusters are detected at least in I-band, only the information 
about individual cluster ages and extinction allows to avoid 
uncertainties of orders of magnitude in star formation rate es- 
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* Based on observations collected at the European Southern 
Observatory, Chile, programme identification numbers 63.N-0528, 
65.N-0577, and 67.B-0504 



timates determined from optical fluxes. From the distribution 
of individual cluster extinction vs. age, which is significantly 
higher for clusters below 8-9 Myr than for older clusters, we 
infer that this is the time by which a typical cluster blows free 
of its native dust cocoon. 
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1. Introduction 



NGC 4038/4039 - "the Antennae" - is one of the closest 
(D ~ 19.3 Mpc, Ho = 75 kms^'Mpc -1 ) examples of merging 
spiral galaxies and thus has been the target of detailed stud- 
ies in many wavelength ranges (e.g . , Whit more et al~l Il999l 
lFabbianoll200ll IWilson et al. ll2Q0lV etc), some aimed at un- 
veilling the star formation properties. However, the term "un- 
veil" is appropriate also in a literal sense, since many of the ac- 
tive star formation sites - especially in the region where the two 
galaxies seem to overlap (therefore often referred to as "over- 
lap region") are hidden by up to several optical magnitudes of 
extinction. Mid- and far infrared wavelengths are essentially 
not affected by this, but the spatial resolution of the images is 
too low to resolve the star formation regions into single star 
clusters. We therefore aimed for a compromise by observing 
the merger in the near infrared, where the extinction is only 
-10% of that at visible wavelengths, but where ground-based 
telescopes currently deliver sub-arcsecond resolution if atmo- 
spheric conditions are favourable. We combine the information 
from these broad band (Ks, ISAAC on the ESO-VLT) and nar- 
row band (CO 2.34jum, ISAAC and Bry, SOFI on the ESO- 
NTT) images with archival Hubble Space Telescope WFPC2 
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Table 1. Summary of the NTT-SOFI and VLT-ISAAC NIR 
imaging observations. 



images (broad-bands roughly corresponding to U, B, V, I and 
narrow band Ha, see lWhitmore et alTl 999: hereafter W99) for 
our investigation of the properties of the numerous compact star 
clusters in the Antennae. 

What we determined from these data are a number of prop- 
erties of the single star clusters, namely their ages, the extinc- 
tion towards them, and their photometric masses. Additionally, 
we derived properties of the ensemble of star clusters as a 
whole: The age distribution over the last ^200 Myrs (includ- 
ing the mass production in the burst), the average extinction 
for cluster populations of different ages, and various other esti- 
mates. The analysis is complemented by an estimate of the bias 
introduced by using only the Ks-band detected clusters, be it an 
age or mass selection effect. 

2. Observations and data reduction 

ISAAC (VLT-ANTU) imaging of NGC 4038/4039 was per- 
formed in ON/OFF mode during the nights 15.04.2001 (Ks- 
band) and 16.04.2001 (CO-band-head filter). The target fit 
completely onto the detector (071484/pixel, total field size 2'5 
x 2'5). Seeing was excellent during both of these photomet- 
ric nights (PSF FWHM of the co-added frames is below 074), 
and the total on-source time was 360s (Ks) and 480s (CO), re- 
spectively. No narrow band continuum observations were per- 
formed, therefore the Ks-image was used for continuum sub- 
traction. 

SOFI (NTT) imaging during 11./12.05.99 covered all NIR 
broad- and several narrow bands, however, seeing was much 
worse than for the ISAAC data (between 076 and 172 PSF 
FWHM), so that only the 077 Bry-data will be used for the 
analysis. Here the field size was roughly twice that of ISAAC 
(4 .'9 x 4.'9, with a pixel size of 07292). Therefore, to save ob- 
serving time, an on-chip offset pattern was observed instead of 
the usual ON/OFF mode. Due to a slight haze, conditions were 
not photometric on either of the nights. 

Table ^ lists descriptions of the observations analyzed in 
this paper for both, the ISAAC and the SOFI runs. 

Reduction of the ISAAC data was performed using the 
IRAF package 1 . It included dark and sky subtraction (either 
using the median of several neighbouring sky images or, where 
this led to residuals, doing pairwise subtraction), and flat field- 
ing by a normalized median of all sky frames. All of the ON 
frames were slightly offset with respect to each other, in order 

1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



to minimize the effect of pixel defects. Therefore, they had to 
be shifted to a common location before using the imcombine 
task (setting the minmax rejection algorithm to reject the high- 
est and the lowest pixel) to combine the single frames. A pho- 
tometric standard (GSPC S279-F, Ks-magnitude 12.031) was 
used for flux calibration (the resulting zero-point was 24.28 
mag). The data reduction procedure for the CO narrow band 
images was identical. 

The resulting Ks-band frame is displayed in Figure^ with 
circles indicating the point-like sources that were detected us- 
ing daofind from the DAOPHOT package. A detection thresh- 
old of 3<r was chosen, leading to 1072 detected objects, which 
are marked in that Figure. 

Reduction of the SOFI data followed essentially the same 
procedure, with the two differences 1) that there were no real 
OFF frames, therefore several images with observing times 
closest to the currently treated frame were medianed for sky 
subtraction. 2) due to a wrong setting of the "Pupil Rotation" 
switch, the object orientation varied from frame to frame, espe- 
cially during the transit of the object. This had to be accounted 
for before the co-addition of the frames. They were rotated by 
the angle derived from the header information. 

HST WFPC2 images were reduced using the IRAF sts- 
das.hst_calib.wfpc package, and photometric calibration used 
the zero-points and colour transformations given in Holtzman 
et al. (1995) in their table 7 and equation 8. For Ho-, we set 
the zero-point to reproduce the magnitudes of several stars and 
clusters lacking Ha emission (which we know from our 3D 
iMeneel et all d200lh and ISAAC spectroscopy iMeneel et al. I 
(2002) data) for a linear interpolation between V and I magni- 
tudes at the central wavelength of the Ha filter (18.316 mag). 

However, in order to be able to combine the optical and the 
near-infrared data sets, we binned the WFPC2 images down 
to the pixel size of the ISAAC frame and rotated them to 
match orientation. Therefore, in the following, whenever the 
unit "pix" is used, it refers to the ISAAC Hawaii array imaging 
scale of07148/pix. 

Photometry was then performed on the target list which 
was created using the Ks-band image. Running daofind on the 
V-band image results in more than 10,000 detected point-like 
sources. However, the analysis of the HST data alone was thor- 
oughly performed by W99. Since we are interested in the ad- 
ditional information brought about by the K-band data, we re- 
strict our analysis to those objects which are detected in Ks. 

In order to estimate the total cluster magnitudes which are 
required for the mass estimates of the clusters, and also in or- 
der to obtain colours combining the optical and near infrared 
bands, we used a curve of growth technique to estimate the 
magnitudes of relatively isolated clusters (distance to nearest 
cluster more than 10 (12) pixels in V (Ks) - band). We assigned 
the magnitude at a radius where the magnitude did not change 
significantly any more. For a number of those clusters we de- 
termined aperture corrections for both bands and applied these 
to the clusters in the more crowded regions, where cluster dis- 
tances were below the mentioned limits. Due to the different 
PSF FWHMs in the two bands (V: 1.6 pix, Ks: 2.6 pix), we 
used different fixed apertures (V: r = 3 pix, Ks: r = 5 pix), 
and the aperture corrections were 0.7 mag and 0.5 mag for V 
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and Ks, respectively. So, in summary, for the isolated clusters 
we used the curve-of-growth technique, and therefore variable 
apertures, to obtain the total magnitudes, whereas for clusters 
with nearby neighbours, we performed aperture photometry 
and applied an average aperture correction. 

We then determined optical colours (U-B, B-V, V-I) all 
within the same apertures (r = 3 pix), since the PSF did not 
vary significantly over the optical wavelength range. 

For the computation of Ha fluxes, we used the magnitudes 
of the Hff photometry within the r = 3 pix aperture and sub- 
tracted a continuum flux estimated from a linear interpolation 



of the V- and I-band fluxes at the central wavelength of the used 
narrow band filter (A cen = 658 nm). We decided to take such a 
small aperture to avoid including flux from neighbouring clus- 
ters in the cases of crowded regions. Radial plots performed 
on some of the isolated, bright clusters in the continuum- 
subtracted He image showed us that most (70%-95%) of the 
emission line flux was indeed contained within this aperture. 
However, in the crowded regions with many clusters around 7 
Myr, overlap in the emission l ine flux is quite c ommo n (see 
the red channel of Figure 4 in Whit more et al. 1 111999)). The 
largest shells of some older clusters have radii of more than 20 
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Fig. 2. V-I colours in the Johnson-Cousins system and the dif- 
ference between the Cousins and the Johnson I-band magni- 
tudes are plott ed for stars of s pectral types B2 t o K4, taken from 
iMoffett & Barnes! dl979l) and Landolt) dl983l) . Over-plotted is 
the best fit to the data, which is represented by the given ana- 
lytic expression. 

pixels, and we clearly cannot integrate out to these distances 
without including a lot of neighbouring clusters. Our selection 
of a 3 pix radius is therefore the attempt to optimize the result 
for young clusters with strong emission. For older clusters with 
extended, low-luminosity shells, this will lead to age estimates 
which are too old, but we will only use the information from 
the Ho- emission for strong Ho- emitters, anyhow. 

Bry-fluxes were more difficult to obtain, due to the much 
lower spatial resolution. Here we proceeded as follows: We 
convolved the Ks-band image with a Gaussian, so that the PSF 
FWHM matched that of the Bry image. Then we performed 
photometry on both images (the zero-point of the Bry image set 
such that the K-band magnitudes of stars and clusters known to 
lack Bry emission 2 were reproduced) for the 1072 point source 
positions detected in the higher resolution image. We then used 
the Ks-band fluxes as the continuum fluxes, and determined 
Bry fluxes and equivalent widths from the difference of the 
two images. In the locations of extremely strong Bry emission 
(equivalent width above 150 A), the Bry flux is likely to be un- 
derestimated by 5-7%, due to the line contribution to the con- 
tinuum filter. 

CO (2.34yum) images had the same PSF FWHM as the K- 
band images, therefore the procedure was facilitated here com- 
pared to the SOFI narrow band data. We set the zero-point of 
the image to reproduce the magnitudes of stars and isolated 
clusters known to lack CO band-head absorption 2 (Z = 21.475 

2 from our 3D (Mengel et al. 1 1200 ll) and ISAAC spectroscopy 
I Men gel et al. 12003) data 
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Fig. 3. Two curves which display the transmission of the CO2.34 
narrow band filter and of the Ks broadband filter installed 
at VLT-ISAAC. The spectrum is a red supergiant of spectra l 
type M3Iab from the stellar library o flLancon & Wood! (2000 ), 
which was shifted in wavelength appropriate for the redshift of 
the Antennae galaxies. Note that the Ks filter cuts off just short 
of the CO band-head and is not substantially influenced by its 
absorption. 
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Fig. 4. The CO index as a function of effective temperature 
estimated for the redshifte d spectra from the stellar library 
of lLancon & Wood 1 12000) using the ESO filter transmission 
curves from Fig. [3] The lines are the two best linear fits, for 
supergiants and giants, respectively. 



mag). The CO index is created taking the ratio of the flux in 
the CO narrow band filter with that in the Ks-band filter and 
expressing it in magnitudes. Using the broad-band filter for an 
estimate of the continuum should be fine, because it cuts off just 
short of the rest wavelength of the first CO band-head (see also 
Fig. [3}- Only in case of extremely high Bry equivalent widths 
(above 150 A), a CO index of up to 0.07 mag can be obtained 
despite the absence of CO absorption. The aperture used for 
determination of the CO index was 5 pix. 
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Fig. 5. The evolution of the CO index using the ESO CO2.34 and 
Ks filters, as it is predicted by Starburst99 after applying cor- 
rections for the definition of the index. The general behaviour 
is similar to other CO indices or equivalent widths definitions 
with only the absolute values differing. 

3. Analysis and Results 

We used several different approaches to estimate the ages and 
extinction in individual clusters based on broad- and narrow- 
band photometry (we are only talking about "clusters", because 
even the brightest supergiants with M K --12 mag are unfor- 
tunately still below our detection limit). Subsequently, these 
methods were combined to yield the cluster properties. Broad- 
band photometry was used to simultaneously estimate cluster 
age and extinction by fittin g the (usually) five dat a points to 
colours from a Starburst99 jLeitherer et al.1ll999l) model for 
ages up to 500 Myrs. 

In order to be able to compare our photometry with the 
predictions of Starburst99, we applied three corrections to the 
models: 

- In V-band and Ks-band, we accounted for line emission in- 
cluded in the filter bandpass. 

- The filter systems on board HST do not exactly match the 
assumed filter transmission in the Starburst99 code, but dif- 
ferences are only significant for I-band. We therefore trans- 
formed the theoretical I-magnitudes to the HST photomet- 
ric system. 

- The CO index computed by Starburst99 is a spectroscopic 
index which does not exactly correspond to the index that 
results from our choice of filters, so we derived a new ana- 
lytic expression for the variation of CO index with effective 
temperature and inserted this in the Starburst99 code. 

The procedures we followed were in detail: 



3.1. V-band and Ks-band: inclusion of line emission 

The HST broadband filter that most closely resembles the 
Johnson V-band filter is F555W, and the zero-points and trans- 
formations described in section |5] actually convert the mag- 
nitudes to the Johnson system. The transmitting wavelength 
range of this filter includes several emission lines, amongst 
them Ha and H/3, 

We made a crude attempt at taking this line emission into 
account in the theoretical predictions of the Starburst99 code. 
The H/3 line will contribute the largest amount of hydrogen re- 
combination line emission, because at the position of the Ha 
line (6563 A), the filter transmission is already down at the 
«10% level, while the intrinsically fainter H/3 line (4861.33 
A) lies in the region of highest transmission (*90%). Also the 
[OIII] line lies in the region of high transmission, and even 
though its strength depends on environmental parameters, it is 
typically much stronger than H/3. We included its contribution 
by using the "typical G alactic HH- region " line ratios listed in 
lAnders and Fritze-v. Alv ensleben (2003), amounting to a fac- 
tor of 5.5. To determine the absolute line fluxes we used the V- 
band magnitudes and the Ha equivalent widths predicted by the 
Starburst99 code, for the strength of the H/3 line we assumed 
the the oretical line ratio o f Ha/H/3 of 3.3 for Case B recombi- 
nation ( Osterb rock II 1 9741) . Modifications in the predictions of 
V-band magnitude are only noticeable for cluster ages younger 
than ^6 Myrs. In a similar way, Bry emission was included 
in the Ks-band filter predictions. However, the contribution is 
quite small. Only for clusters below 4 Myrs, the difference in 
Ks-magnitude amounts to more than 0.05 mag. 

3.2. I-band photometry 

While the HST filter set most closely resembles the Johnson- 
Cousins-system, which means that the I-band filter is in the 
Cousins system, the Starburst99 code assumed an I-band fil- 
ter transmission curve corresponding to the Johnson system. 
In order to obtain a transformation between the two systems, 
we used photometry in both systems f or stars of spectral types 
B2 to K 4 taken fromlLandoltl dl983h for the Cousins system 
and from Moffett & Barnes ( 1979) for the Johnson system. We 
converted the Starburst99 I-band magnitudes from the Johnson 
system to the Cousins system using the following expression: 

^Cousins ~ ^Johnson 0.28 * (V /) Johnson-Cousins ~ 0.01 

Figure |3 shows the data points we used together with this 
best fit. 

3.3. CO index 

As mentioned in section[2] we used the narrow band C02.34//m 
filter and the broadband Ks filter for estimating the CO in- 
dex. This arrangement is a c lose approximation to the pho- 
tometric CO index defined in Klein mann & Half! {T986) who 
use another narrow band filter for an estimate of the K- 
band continuum. The spectroscopic CO index which is com- 
puted in Starburst99 is the spectroscopic CO index defined by 
Origlia & Oliva I |2£)()0). Rather than attempting a conversion, 
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Fig. 6. Spatial distribution of the three narrow band age diagnostics, with the colour coding as indicated in the individual plots. 
Generally, blue and cyan colours indicate younger clusters than magenta and red, but the legends for each plot should be con- 
sulted for individual comparisons. The general distributions of ages are very similar, and the ages determined from the various 
diagnostics coincide within a few Myr for most clusters. There are a few exceptions, however, where one or more of the indicators 
show a deviant value. 



we used the giant and superg iant stars in the near infrared stel- 
lar library assembled by Lan con & Wo od (2000). The K-band 
spectra of the stars were redshifted approximately to the red- 
shift of the Antennae galaxies (z«0.005) and we estimated their 
CO indices using the ESO filter transmission curves for the Ks 
and the CO2.34 filter. From this analysis, we determined two 
analytic expressions for the CO index as a function of effective 
temperature, one for the giant and one for the supergiant stars. 



We did not determine the CO indices for dwarf stars, because at 
cluster ages below a few hundred million years, their contribu- 
tion is negligible. We note that the library of lLancon & Woodl 
(2000) only included high luminosity stars, but had we ex- 
pected a substantial contribution for our purposes, we would 
have included them using a different library. 
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The two analytic expressions we derived are: 



COindex = 0.645 - 0.942 ■■ 



COindex = 1.008 -2.186 = 



'eff 



10000£ 



'eff 



10000K 



(giants) 



(supergiants) 



and the data points and the corresponding fits are shown in 
Fig. El 

W e replaced the corresponding expressions (IDovon et al. I 
1 19941) in the Starburst99 code by these two expressions, then 
reran it to produce the evolution of the feature for an instanta- 
neous burst at solar metallicity for an IMF with Salpeter slope 
between 1 and 100 M©. The result is shown in Fig. [5] The gen- 
eral evolution of the feature is the same as for other age tracers 
using this spectral feature - like the above mentioned photomet- 
ric and spectroscopic indices, or the equivalent width. Only the 
absolute values have changed, and they can now be compared 
with our observations. See section l3~5l for this comparison. 

For six clusters, we have information about the equivalent 
width of the CO band-head (as defined in Or iglia et al. II 19981 
and used in Starburst 99) from spectroscopic or integral field 
spectroscopy data ("see lMengel et al. I200ll l2002). We compare 
(see Table |2j these values and the corresponding age predic- 
tions to the CO index and assigned ages for two reasons: Firstly, 
in order to check the reliability of the technique, and secondly, 
to assess how strong the influence of extinction on the observed 
CO index is. We ended up extinction-correcting the CO index 
according to the following formula, which makes a noticeable 
difference only for clusters with an Ay of above »2 mag. 

COindex(extinctioncorrected) — COindex + 0.012 * Ay 

This accounts for the differential extinction expected from 
the difference in central wavelength between the two filters. 
The two higher extinction clusters (W95-355 and W95-80) in 
Table |3 would allow for an even higher correction factor, but 
we stuck to this physically justified correction. 

As can be seen from Table [2] the agreement is quite good 
for these clusters, but it is also obvious from the two clusters 
with high extinction that it improves agreement to apply an ex- 
tinction correction for Ay above * 2 mag. It is unfortunate for 
this purpose that our spectroscopic studies were heavily biased 
towards clusters around 10 Myrs of age, because this narrows 
down the age range available for comparison between spec- 
troscopy and narrow band imaging. 

3.4. Our Standard Model 

We used the Starburst99 dLeitherer et al. 1 1 19991) model, be- 
cause it focuses on young stellar populations, which are dom- 
inant in this young merger. Nevertheless, we would have liked 
a model which also takes into account the thermally pulsing 
AGB phase, to be able to assess the impact they have on K-band 
luminosi ty. It is large according to some authors (see for exam- 
ple Forster-Schreiber et al. 2003; and references therein, and 
Schul z et al. N2002l) 1. but depends a lot on assumptions like the 
duration of that phase. Just before this pa per was ready for sub- 
mission, the new version of Starburst99 dVazauez & Leithererl 
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2005; including the Padova tracks with full AGB evolution) be- 
came available. Too late to re-do the analysis, but in time to add 
one track to Fig. EH It shows that the impact on our interpreta- 
tion would have been marginal if we had used the new models 
instead. The 10 Myr (red supergiant) peak in K-band luminos- 
ity has a smaller amplitude and is somewhat wider in the new 
model (which would probably have caused our age distribution 
to get slightly wider). We emphasize however that the peak in- 
troduced by the population of TP- AGB stars of around 200 Myr 
is not pronounced enough to have shifted a significant fraction 
of the K-band selected clusters into that age range within our 
detection limit. Thus the differences in the models did not lead 
to a significant number of incorrect (younger) age assignments. 

The characteristics of the model we used throughout this 
paper (unless noted otherwise) was an instantaneous burst, with 
solar metallicity, and a Salpeter IMF with a low and high mass 
cut-off of 1 and 100 M© respectively. These parameters are jus- 
tifiable. The instantaneous burst is applicable because the star 
clusters are typically very small (effective radii of a few par- 
sec) and very concentrated, which means that the star formation 
event is essentially instantaneous. 

Solar metallicity or just slightly above is what we deter- 
mined for a few bright individual clusters in the Antennae from 
an analysis of UVES high spect ral resolution spectra (paper 
in preparation). Note that Vazauez & Leitherer (2005) caution 
against using their model (or any current evolutionary synthesis 
model) for stellar populations with sub-solar abundances and a 
significant population of red supergiants. This is the main rea- 
son why we decided to fix the metallicity at a reasonable av- 
erage value, rather than making it a fit parameter, d espite the 
necessary caution suggested bv lAnders et al. I ll2004l) . 
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Identification 


A v 


Wco 


CO index (ext.corr) 


Age (Wro) 


Age (CO index) 




Imael 


rAi 


Tmaffl 


rMvrsl 


TMvrsl 


W99-1 


0.094 


17.5 + 1 


0.199 


8.6.. 12.7 


8.6.. 12.5 


W99-2 


0.03 


16.2 ± 0.2 


0.165 


8.5 or 12.8 


8.4 or 12.7 


W99-15 


1.02 


17.0 ± 0.2 


0.279 


8.6.. 11.5 


8.7.. 12.3 


W99-16 


«o 


19 ±4 


0.673' 


8.2.. 12.9 


8.7.. 12.3 


WS95-355 


2.88 


16.3+0.2 


0.134 


8.5 or 12.8 


8.2 or 13 


WS95-80 


4.01 





-0.02 


<7.3 


< 7.3 



Table 2. Comparison of ages derived from Wco obtained from spectroscopy (3D at the AAT and ISAAC at the VLT) with our 
narrow band determined CO index. Agreement is generally quite g ood. For definition of Wr r>, C O index and extinction correction , 
please see text. The identification refers to the numbers used in IWhitmore et al. I J 1994 ' and Whitm ore & Schweizer I |l995), 
respectively. 1 This value is completely off-scale, even if uncertainties are large. The age we assign is the range corresponding to 
the maximum CO index of *0.22. 



The IMF is pro bably only on average around Salpeter 
(Men gel etal1 bo02). and we expect it to vary for the whole 
cluster population in a si milar manner to wha t was seen for 
the subsample analyzed in iMengel et al.¥ 2002). Since we have 
not yet found a parameter which to correlates significantly with 
IMF slope or l ower mass cutoff, we assume this average value 
determined bv IMengel et al. I d2Q02l) . Note that assuming, for 
example, a lower m ass cutoff of 0. 1 Mq for the Salpeter IMF, 
or a Kroupa IMF (Kroupa 2001) instead, does not change 
our conclusions derived from relative quantities like colours or 
equivalent widths, because those IMFs largely add mass in the 
form of low-mass stars, which do not contribute a significant 
amount of light at young ages. Only the cluster mass derived 
from the photometry would be higher than the ones we obtain, 
therefore our masses represent lower limits. For example, the 
mass of a cluster for a Salpeter IMF between 0.1 and 100 M0 
is higher by a factor 2.6. 

3.5. Age Determination 

Our data set includes a number of age indicators, therefore we 
used the following strategy to assign individual ages to the Ks- 
band detected clusters. We deemed this necessary (rather than 
a simple averaging of all available ages) because of the on/off 
nature of the narrow band diagnostics: For example, absence 
of CO absorption in a sufficiently bright cluster constrains the 
age to be anything below 6.3 Myrs, which can be used to verify 
ages from other diagnostics, but not for straight averaging. In 
essence, what we did was: 

- Firstly, we have the broadband data covering U, B, V, I, 
and Ks, and in principle suitable to age-date clusters of all 
ages (the large wavelength coverage, especially the inclu- 
sion of Ks-band photometry, is suitab le for break i ng the 
age-extinction-degeneracy - see also And ers et al. I ([2004) 
for a more general discussion). 

- Secondly, we can use narrow band emission line images: 
Ha and Bry. The presence of hydrogen recombination line 
emission indicates a very young age of the cluster, below 

Myrs, and we use the equivalent width to determine the 
age if it is below this limit. 



- Thirdly, CO band-head absorption narrow band images, 
where the detected spectral feature reveals the presence of 
red supergiants or giants, and is most efficient in identifying 
star clusters around 10 Myrs of age. 

As mentioned, we used the broadband data to break the de- 
generacy in age and extinction. To do that, we performed a si- 
multaneous fit of the theoretical spectra for ages up to 500 Myrs 
to the observed broadband magnitudes. The amount of extinc- 
tion was determined by reddening the theoretical spectrum ac- 
cording to the Mathis (1990) extinction law. Note that we are 
assuming a foreground screen extinction. We believe this to be 
appropriate, because after the onset of stellar winds and super- 
novae, remnant gas and dust are expected to b e driven out o f the 
cluste r within less than a crossing time (see B oilv & Kroupa I 
120031 which for the clusters we have studied is 5510 5 yr. This 
material forms a cocoon around the cluster which can be de- 
scribed as a foreground screen. A mixed model would have to 
be assumed if we analyzed the average extinction in the galaxy 
from integrated emission, but not in determining the extinction 
towards individual clusters. To simultaneously determine age 
and extinction, we used a simplex algorithm to tour data space, 
with ax 2 minimization technique to find the best fitting pair of 
extinction/age. 

Ages and extinction so obtained were complemented for 
the brighter clusters (only brighter clusters because the narrow 
band images were shallower) by extinction estimates derived 
from the ratio of Ha/Bry flux compared to the value expected 
for a Case B recombination scenario, and by age estimates from 
the Bry equivalent width , the Ha equivalent width and the CO 
index. 

The average ages were usually assigned by averaging the 
available data, but we kept track of the number of available data 
points and their level of agreement by assigning an "age qual- 
ity" parameter, which is highest if all four indicators (broad- 
band fit, W//q,, Wflry, CO index) are present and agree, and is 
lower for fewer data points or disagreement: 

- If all four indicators yielded ages (74 cases), the common 
age was the average of the four, if the age difference be- 
tween any of them was smaller than 4 Mys. The "age qual- 
ity" was set to 4. If only three indicators agreed, their aver- 
age was used and the fourth value ignored, with the quality 
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Clusters with discrepant ages in different indicators 



Age < 6 Myrs 
6 Myrs S Age < 9 Myrs 
9 Myrs g Age < 13 Myrs 
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hi 



6 6.5 7 7.5 

log (Age [yrs]) 



Fig. 8. Colour (age) and size (luminosity) distribution of clus- 
ters that had discrepant age determinations. Roughly 180 clus- 
ters fall into this category. The colour coding represents the 
best age fit determined from using only broad-band colours. 
A 2-dimensional Kolmogorov-Smirnov test indicates that the 
spatial distribution of these clusters is marginally different from 
the distribution of all clusters (probability of 4.6%). This is be- 
cause the clusters that lie in regions of high background relative 
to their total brightnesses are more likely to have discrepant 
ages than the general population of clusters. Thus bright bright 
clusters in the nuclei, for example, tend to have inconsistent 
age determinations, whereas faint clusters with discrepant age 
determinations tend to be distributed more like the general pop- 
ulation of clusters. This hypothesis is confirmed by running 2- 
dimensional K-S tests separately on the bright and faint sub- 
populations with discrepant ages compared to the general pop- 
ulation of bright and faint clusters respectively. The dividing 
magnitude was chosen to be m^=15.5. 



6.5 7 7.5 

log (Age [yrs]) 



value set to 3. If two or more ages were discrepant with re- 
spect to the others, the common age was set to the average 
of the Ha- and the Bry-age, but the quality-value was set to 
a low value of 2. 
- Similar strategies were followed for average age determi- 
nations if fewer than four indicators were available (three: 
376, two: 320, one: 262), essentially averaging the avail- 
able ages and assigning data quality parameters depending 
on the level of agreement. What should be mentioned is that 
for the cases where only one age was available (that would 
usually be the value from the broadband fit), that age was 
also assumed for the average age, and that cases where only 
one or two indicators were available, but where the other 
indicators are expected to be missing were also assigned 
high age quality parameters. For example, a cluster with an 
age determined from the broadband fit and the CO index to 



Fig. 9. Histogram of numbers of star clusters observed in each 
age bin. All clusters are young, and most clusters are either 
around 10 Myrs or around 6 Myrs old. Different colours/line- 
types differentiate between V-selected brightness limits (top) 
and Ks-selected (bottom). The strong peak of clusters of age 
10 Myrs is a selection effect introduced by detecting clusters 
in Ks-band, which has a strong peak in absolute luminosity 
around that age, caused by high-mass stars entering the red su- 
pergiant phase. This means that many low-mass clusters are de- 
tected in that age bin which would be below the detection limit 
at other ages (see also Fig. 1101 . Note the increasing relative 
number of young clusters («4 Myrs) when going to brighter 
clusters in both figures, which is already an indication that more 
clusters are present at younger ages, an even more convincing 
evidence for this being displayed in Fig. II II 
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be older than 10 Myrs is not expected to show Ha or Bry. 
Therefore, it is assigned an age quality parameter of 3.9. 

Out of our 1072 Ks-band detected clusters, 40 could not be 
assigned any ages, 997 had ages from the broadband fit, 590 
from the CO index, 128 from Ws r y and 494 from Wh«- In gen- 
eral, the agreement between the age indicators was good: ex- 
cellent agreement was seen for « 340 clusters, another w 500 
showed good agreement, and only as 230 were vastly discrepant 
(which includes the 40 clusters without an age determination). 
Their distribution is shown in Fig. [8] Most of the brighter ones 
of these objects with bad age determinations were either fore- 
ground stars, or located in the nuclear regions, which suggests 
that for the latter the high and variable background of the nu- 
cleus might be responsible for photometric errors that led to 
wrong age assignments - for example because line emission 
was erroneously detected, or subtracted from the source be- 
cause it is too extended. The fainter objects were distributed 
throughout the merger, which means that a disagreement in age 
determinations for them is probably caused by statistical uncer- 
tainties in the photometry. For the following analyses, clusters 
with an "age quality" parameter below 3.0 were not taken into 
account (which means that foreground stars with their bad age 
fit are automatically excluded from the analysis). 

The age distribution (see Fig. [7] and Fig. |5Jl shows two 
peaks, one at approximately 5-6 Myrs, and one around 10 
Myrs. Older clusters are rare, there are two more peaks at 
around 28 and 200 Myrs, with roughly 20 members each. The 
strong peak around 10 Myrs is a selection effect which arises 
from the very high total Ks-band luminosity around that age, 
leading to detection of large numbers of low mass (between 
10 3 and a few times 10 4 Mq) clusters which are below the de- 
tection limit at other ages. 

However, one interesting aspect of the age distribution is 
that the relative numbers of 6 vs. 10 Myr old clusters shifts with 
limiting magnitude: Going to brighter clusters, the younger 
ones become more dominant. This is expected from the lumi- 
nosity evolution in V-band, where th e youngest clust e rs ar e the 
brightest (see Figs. 47b and 51b in Leithereret al. 1 ll 19991) or 
Fig El m this work for the luminosity evolution in V and K), 
but the absolute Ks-band magnitude peaks at a cluster age of 
«10 Myrs, where the luminosity contribution from the super- 
giants dominates that band. Therefore, we expected to see this 
trend in the plot of the V-band selected clusters (Fig. [5] top), 
but were surprised to see this, even though to a lesser extent, 
also for the Ks-band selected clusters (Fig. [9j bottom). Since 
we used the extinction corrected magnitudes, it is not extinc- 
tion which causes this effect. It rather means that actually more 
high luminosity (and therefore high mass) clusters are present 
at ages around 6 Myrs, compared to those around 10 Myrs. 
This conclusion is also supported by Fig. \^\ which clearly 
shows the drop in cluster numbers between 4 and 15 Myrs 
for the high mass clusters. We performed a one-dimensional 
Kolmogorov-Smirnov (K-S) test in order to assess the signifi- 
cance of the decrease in cluster numbers. For the clusters above 
1O 5 M0, we created many hundreds of random comparison pop- 
ulations, one assuming a constant probability in time, the other 
assuming an exponential probability distribution, both with the 
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Fig. 10. Two different representations showing the complete- 
ness and/or detection limits for V and Ks band for our clus- 
ters. The upper graph shows completion limits for clusters with 
different total masses in V- and Ks-band (both lower curves 
are for a lxlO 5 M0 cluster, for Ks-band, we added 5xl0 5 M© 
and lxl0 6 M Q ). This shows that in V-band, 50% of all clus- 
ters above a mass of lxl0 5 M o can be detected up to an age 
above 1 Gyr, whereas this cluster would drop below the Ks- 
band completion limit already at 25 Myrs of age. The red 
(dot-dashed ) line was, like the o ther graphs, created using a 
Starburst99 JLeitherer et al. 1 1999) model (see text for our stan- 
dard assumpt ions and modi fications), but this time for v5.0 
dVazauez & Leitherer 1120051) . which, amongst other modifica- 
tions, includes the TP-AGB phase, which is responsible for the 
bump after 10 8 years. The lower graph is in support of our 
claim that many more clusters exist with ages below 25 Myrs, 
and that this is not a simple selection effect, because it shows 
that clusters above a certain mass (but for comparison needs to 
be above the completeness- or detection limit for all ages then) 
are much more numerous below 25 Mvrs than above that ase. 
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Fig. 11. This plot shows the distribution of clusters above cer- 
tain mass limits (10 5 and 10 6 Mq, respectively), as a function 
of age up to an age of 25 Myrs, which is the age where a clus- 
ter of 10 5 M0 reaches the Ks-band completeness limit. This 
means that in our plot, essentially all clusters above the given 
mass limits should be detected, if they are present. The strong 
decline in cluster numbers up to an age of 15 Myrs therefore 
means that indeed many more clusters are present around 6 
Myrs, meaning that either more clusters were produced or that a 
large fraction has been destroyed on timescales of 10-15 Myrs. 
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Fig. 13. Spatial distribution of extinction values (Ay) deter- 
mined for the clusters. Values from the broadband fit and from 
the line ratios were combined to a common average where this 
was possible. Note the excellent agreement between high ex- 
tinction values (red and magenta dots) and the lanes of dust 
obvious in the multi-colour image (the same as in Fig. [21 only 
with different cut levels). 



same total number of clusters as observed. In neither case did 
we concretely account for a selection function because as we 
have shown earlier, our selection should be robust for clusters 
with masses above 10 5 Mo and ages <30-100 Myrs. The prob- 
ability that the observed population was drawn from a parent 
population with a constant distribution of cluster numbers in 
time is virtually zero (only 3.5%). The population is not consis- 
tent with being a simple exponential either - only a probability 
of 56%. For the clusters above 1O 5 M0, the only simple robust 
conclusion we can draw without conjuring up more complex 
models, is that the distribution of cluster ages is consistent with 
our visual impression which is the number of clusters is declin- 
ing significantly with age. 

3.6. Extinction 

A similar procedure like for the ages was followed for the 
determination of extinction Ay (we are assuming foreground 
screen extinction, see section l3~5Y For most clusters, it was 
determined from broadband colours, using the fitting routine 
described in section l3~5l for a Mathis (1990) extinction law, but 
for those with m^ < 17.3 mag and Wb, t > 80A, the value 
determined from the observed flux ratios of He and Bry was 



averaged with the broadband determined value by averaging 
the fluxes and converting back to magnitudes. 

As could be expected from just looking at the false colour 
image (TFig. 1121. extinction is very variable within the merger, 
peaking at extreme values of Ay above 10 mag in the over- 
lap region. All extinction values quoted in this article assume a 
foreground screen extinction model. 

The distribution of extinction is shown in Fig.^] where we 
plotted a colour-coded map of Ay determined from the broad- 
band fit, and, where possible, in combination with the line ratio 
value. Underlying t his distribution is the CO map taken from 
IWilsonetal.l fcoOl). to demonstrate the excellent agreement 
between high extinction clusters and density peaks in the CO 
map. 

The average extinction is around Ay = 1.3 mag, and it is 
interesting to analyze the age dependence of the extinction dis- 
tribution, reproduced in Figure [21 for information about the 
clearing time of the clusters. 

The typical extinction decreases with increasing age, as can 
be seen in Fig. [21 which shows the median extinction values as 
a function of age in bins of 2 Myrs, together with the individual 
data points which were used in the median. See Sec.0]for an 
analysis of this evolution. 

These age-averaged numbers artificially increase the ex- 
tinction values for older clusters, at least if the purpose of the 
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Fig. 12. False colour image of NGC 4038/39 with U+B in the blue, V+I in the green (these four images were obtained from 
the HST archive, see W99 for details) and Ks (ISAAC, VLT) in the red channel. Most of the reddening seen in the clusters 
in the "overlap region" between the two nuclei (and elsewhere in the merger) is caused by dust absorption and therefore gives 
a first-glance impression of the patchy and partially high extinction encountered in the merger, which was the main driver for 
obtaining our near-infrared observations. FWHM of the PSF is (X'4 (HST images were smoothed), the region displayed covers 
approximately 2x2 arcminutes, north is up, east to the left. 



analysis is to determine the extinction directly related to the 
cluster in question, for example arising from its own dust co- 
coon. This is suggested by Fig.^J Visual inspection of the lo- 
cations of high- and low-extinction clusters above 9 Myrs with 
respect to the locations of young, high-extinction clusters show 



that the high extinction older clusters appear almost exclusively 
in the direct vicinity of younger clusters, which we took as an 
indication that they only happen to lie within or behind the 
molecular cloud related to the formation of the younger, neigh- 
bouring cluster. 
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Fig. 14. Small black dots are individual Ay, age data points for 
clusters with good age estimates. The large red dots are the me- 
dians of extinction values found in 2 Myr wide bins (the last bin 
shown contains 36 data points, many of which have Ay close to 
and therefore overlap). Despite the large scatter in Ay values 
in any given age bin, the expected trend of lower Ay with in- 
creasing age is obvious, and many of the clusters older than 8-9 
Myrs have very low extinction. This lends further support that 
this is the age where many of the clusters emerge from their 
natal dust cocoon. 

In an attempt to quantify the robustness of this visual 
impression, we calculated (via a two-dimensional K-S-test) 
the probability that the population of clusters that are young 
(age<9 Myrs) and highly extinguished (Ay >1.5mag), old 
(age>9 Myrs) and highly extinguished (Ay >2.0mag), or old 
with low extinction (Ay <0.8mag) are spatially distinct. Our 
visual impression is confirmed. We find that the probability that 
the young and highly extinguished and the old and highly ex- 
tinguished to be spatially distinct to be very low (only 22%). 
Thus it is likely that these clusters are drawn from the same 
sample. However this is not the case for the young and highly 
extinguished clusters and the old clusters with low extinction. 
These samples are have a probability of <<1% of being drawn 
from the same population spatially. 

We plotted Fig. [21 with lower age limits for the "older" 
clusters increasing from 6 to 12 Myrs in steps of 1 Myr (but we 
do not show these plots here), and the described effect becomes 
obvious between 8 and 9 Myrs, which we interpret to be the 
time when most star clusters have cleared away their natal dust 
cocoon. This is the same age range which is suggested by Fig. 

m 

Mostly caused by extinction, a fraction of the Ks-band de- 
tected clusters is not detected at shorter wavelengths. To esti- 
mate the impact extinction has on cluster detection and esti- 
mates of, for example, star formation rates based on observa- 



Fig. 15. Extinction values colour coded in their distribution in 
the Antennae. The 10 Myr clusters with high extinction are 
preferentially located near younger clusters, while those with 
low extinction values are often either single or in groups of 
older clusters. 



tions in certain optical bands, we determined which fractions 
of clusters were undetected, how much mass they contribute, 
and how this compares to the mass which would be derived for 
the total cluster mass if only the observed light in clusters for 
the given band is used for mass determination. 

Table 13.61 shows that the total fraction of clusters unde- 
tected in any given optical band, compared to the Ks-band, is 
rather low, despite the sometimes fairly high extinction, thanks 
to the fact that the HST observations are up to a few magnitudes 
deeper: Even in U-band, only 16% of the clusters are not de- 
tected at all, and those clusters contribute only 12% of the total 
mass in clusters. These numbers become even more favourable 
going to V or I band. 

The impact of extinction becomes obvious, however, if one 
compares the masses that would be derived from optical clus- 
ter light without applying an extinction correction. Using the 
individual cluster ages to derive individual cluster masses from 
just the observed magnitudes, like we did for Ks-band, U and B 
band lead to a total mass estimate which is just a quarter of the 
mass derived from Ks, and this improves only slightly towards 
V and I band, where a third of the Ks-band mass is derived. 

Not taking into account the individual age information, but 
rather assuming average ages and extinction can lead to more or 
less arbitrary results: For example, assuming an average age of 
5 Myrs for the clusters and converting the observed U-band flux 
to a total mass yields 33% of the mass derived from the extinc- 
tion corrected Ks-magnitudes (Ks-mass), whereas assuming an 
average age of 25 Myrs and an average extinction of Ay of 1 
mag gives 9 times the Ks-mass. 
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Band 


number of clusters detected 


mass of clusters not 


mass derived from m Ba „d and 




in Ks, but not in Band 


detected in Band over Ks mass 


age over Ks-mass 


U 


16% 


12% 


26% 


B 


13% 


10% 


27% 


V 


6% 


5% 


33% 


I 


9% 


6% 


33% 



Table 3. This table compares the photometry of detected and undetected clusters in the optical bands to the Ks band, in order to 
assess the impact of extinction on e.g. cluster detection and the star formation rates derived from them. The first column gives 
the names of the bands (rather, the Johnson equivalents to the HST broadband filters), the second column gives the number of 
Ks-detected clusters which are not detected in the specified band (in percent). The Ks-mass, which is referred to in columns 
3 and 4, was obtained in the following way: The extinction corrected Ks-magnitudes for individual clusters were compared to 
the magnitude expected for a cluster of 1x10 6 Mq (using Starburst99 for an instantaneous burst, solar metallicity, Salpeter IMF 
between 1 and 100 Mq). These masses were summed up to derive the total Ks-band mass. Column 3 is the total Ks mass of the 
clusters not detected in the specified band, whereas column 4 specifies the mass of all the clusters detected in the given band, if 
only the not extinction corrected flux in the given band is used (together with cluster age) to determine the mass. 



3. 7. Photometric masses 

Knowing the ages of the clusters is the presupposition for de- 
termining their photometric masses, M p h t- The extinction cor- 
rected Ks-band magnitudes are compared to those expected for 
a cluster of a given age and mass for the model parameters. We 
use Ks-band magnitudes, rather than any of the optical colours 
for two reasons: Firstly, in Ks-band all of our clusters are de- 
tected, so it is not necessary to extrapolate to one of the optical 
colours for an assumed age and extinction, which always intro- 
duces additional uncertainties. And secondly, Ks-band is much 
less affected by extinction, which additionally minimizes the 
uncertainties. 

The Starburst99 JLeitherer et al. l [l999) model parameters 
we used were again instantaneous burst, solar metallicity and 
an IMF with Salpeter slope between 1 and 100 Mq (note that 
for an IMF extending between 0. 1 and 100 Mq, the photomet- 
ric masses would increase by a factor of 2.6). From the result- 
ing individual cluster masses we determine the total mass of 
stars produced in clusters during this starburst and the corre- 
sponding star formation rate (or rather, lower limits to both, 
since there are still many clusters expected to be below the de- 
tection limit), and a cluster mass function which we compare 
to the luminosity function. 

The total mass produced in Ks-band detected clusters with 
an age determination is ^430 x10 6 Mq. The error introduced 
by not taking into account those clusters without age determi- 
nation is fairly small, since it is mostly the faint, low mass clus- 
ters which lack an age estimate. And, as mentioned above, the 
SFR we determine here will be a lower limit. From Fig. [9] we 
see that the bulk of star clusters are younger than 25 Myrs, 
which yields a SFR of SFR c / t(s , e „ * 16 Mq yr _1 , which is al- 
most 20 times the Milky Way SFR, and a moderate to normal 
value for a starbursting galaxy. Interestingly, the total molecu- 
lar gas mass estimated for the Antennae is approximately 10 10 
Mq from CO observations (Wils on et al. 1120071) which is sub- 
stantially greater than the estimated total cluster mass - approx- 
imately a factor of 10-20 (depending on whether the IMF is 
Salpeter from 1-100 Mq or a reasonable extrapolation like a 
Kroupa IMF with a minimum mass of 0.1 Mq) of the cluster 
mass reported here. 



Thus the current gas depletion time would be about 250 to 
500 Myrs, if a star formation efficiency of 100% was assumed. 
For a more reasonable assumption on the star formation effi- 
ciency, like 30%, the depletion time needs to be scaled down 
accordingly, but in any case, this supports further intense star- 
formation over the expected merger time scale. 

Taking into account our conclusions concerning cluster dis- 
solution, and instead assuming constant star formation over the 
last 25 Myr, our lower limit concerning the star formation rate 
increases by at least a factor 2.5, and would then be above 40 
Mq yr- 1 . 

3.8. Cluster Mass Function 

The cluster mass function in merging galaxies is of special 
importance with respect to star cluster formation and evolu- 
tion, because globular cluster mass functions (and since they 
have a small range of ages, also their luminosity functions) 
are Gaussian (or similarly single-peaked), with a typical mass 
around 2x10 5 Mq. However, luminosity functions of clusters 
in merging galaxies show a power law with a typical slope 
a = -2. Cluster mass functions are difficult to determine, be- 
cause, like in this work, it requires the determination of individ- 
ual ages for these young clusters where the intrinsic luminosity 
varies strongly with age. Nevertheless, the mass function for 
the Ant ennae clusters was co nstructed from HST data by Zhang 
& Fall ( Z hang & Fall 1 19991) . and found to be a power law with 
approximately a = -2, with no indication of a turnover at any 
mass. A turnover of power laws is consistent with th e luminos- 
ity fun ction constructed from the same data by Whit more et al. I 
( 1999), with the turnover roughly at 1x10 5 Mq, comparable to 
the characteristic mass of globular clusters mentioned above. 

From our Ks-band data, we created the luminosity- and 
mass functions, and even though the data are shallower in limit- 
ing magnitude, we hope to benefit from the decreased influence 
of extinction on our data. 

Luminosity and mass functions are displayed in Fig.ll6land 

EU 

The extinction corrected Ks-band luminosity function of 
star clusters in NGC 4038/4039 can be fit by two power laws 
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Fig. 16. Extinction corrected Ks-band luminosity function of 
star clusters in NGC 4038/4039. Completeness limit (50% of 
clusters detected) and completeness corrected data points are 
indicated. The data can be fit by two power laws with a turnover 
magnitude of m^ = 15.45. Comparing the corresponding ab- 
solute Ks magnitude of M Ks - -15.96 mag to that of a 1O 6 M0 
cluster of the median age of 9.2 Myrs (Mjf, = -17.52 mag) 
yields a turnover mass of 2.4xlO 5 M0. 



("broken power law") with a turnover magnitude of m Ks - 
15.45, and the slopes a = -2.6 and a = -1.37, respec- 
tivel y. The slopes compar e fairly well with those determined 
bv lWhitmore et al"l i 19991) from the objects located on the PC 
chip which had from an analysis of their Q-parameters been 
confirmed to be star clusters. This fact is somewhat surpris- 
ing, because the turnover masses in the Ks- and V-luminosity 
functions need not necessarily be the same, since the cluster 
luminosity function is a function of cluster ages (and therefore 
mass-to-light-ratios) and masses. 

Comparing the absolute turnover Ks magnitude of M/j s = 
-15.96 mag to that of a 1O 6 M0 cluster of the median age 
of 9.2 Myrs (Mk s = -17.52 mag) yields a turnover mass of 
2.4xlO 5 M0. This could be interpreted as the progenitor to the 
typical mass peak in the globular cluster mass function. 

In Fig.[nJ the cluster mass function can also be fit by a bro- 
ken power law, with the two slopes of a = -1.9 and a = —1.2, 
The turnover mass here is around 3xlO 5 M0, which is in rea- 
sonable agreement with the value derived from the luminosity 
function. 

However, as we will discuss in Sec. |4] we do not want to 
over-interpret this result, since both random and systematic un- 
certainties are rather large, and the latter difficult to overcome. 
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Fig. 17. Mass function of Ks-band detected star clusters in 
NGC 4038/4039, and fits using two power laws (slope -1.9 and 
-1.2). Masses were determined using individual cluster ages 
and comparing the observed magnitudes to those predicted by 
Starburst99 for a 1O 6 M0 cluster. The completeness limit is not 
exact, because it is the completeness limit in magnitudes (see 
Fig. I16> which was converted to a mass for the median clus- 
ter age of 9.2 Myrs. For some ages which are quite well rep- 
resented, for example around 6 Myrs, the completeness limit 
lies rather around 1O 5 M0, which makes it really difficult to ap- 
ply a completeness correction for this plot, which we there- 
fore omitted. Also these data can be fit by two power laws, 
with a turnover mass around 3xlO 5 M0, in reasonable agree- 
ment with the number derived from the luminosity function. 
The assumption that a large fraction of the highest mass clus- 
ters could be composed of cluster complexes, rather than indi- 
vidual clusters, could be excluded by visual inspection of both, 
the K-band image and the higher resolution I-band image at the 
locations of the high-mass clusters. At least down to the spatial 
scales which can be resolved, less than 5% of the clusters with 
masses above 1O 5 M0 look like they might be affected by mul- 
tiplicity. However, very close binary clusters (like for example 
NGC 1569 A, which could be resolved into two com ponents 
with a separation of 072 on HST images ("see lde Marchi et al"l 
1997)) cannot be excluded for the Antennae clusters, which lie 
at almost 10 times the distance of NGC 1569. 

3.9. Bolometric luminosity 

We used the combination of individual photometric masses 
and ages to determine the bolometric luminosity expected to 
be emitted from a given cluster, using the predictions from 
Starburst99 Leitherer et al. where the bolometric lu- 

minosity is obtained by integrating the spectral energy distri- 
bution (without the nebular continuum). From this, we created 
an artificial "bolo metric image" to compare it to the ISOCAM 
15/mi image (Mirabel et al. 1998; see also the comparison of 
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Fig. 1 8. Overlay of ISOCAM 15yum image jMirabel et all 
1998) with simulated mi, i contours. The simulated mt, i im- 
age was created by calculating the expected bolometric mag- 
nitude of each Ks-detected cluster and feeding this, together 
with the position, into the IRAF daophot.addstar routine. The 
resulting image was then Gaussian smoothed to have a resolu- 
tion comparable to the ISOCAM 15/mi image. In general, the 
good agreement between the two maps is noticeable, but apart 
from that, we note that the brightest 15/mi cluster (indicated 
by the arrow at the bottom) does not have the highest bolo- 
metric luminosity, and we identify two other regions where the 
simulated bolometric luminosity is substantially lower than ex- 
pected, compared to the 15//m flux. Those three regions are 
indicated with the arrows. 



15pm image and CO map in Wi lson et al. 2001), t o the radio 
continuum emission observed bvlNeff & Ulvestad Id2000t) . and 
to CO emission dWilson et al.l200ll) . 

Our main goal here was to see whether the mass and age 
of the very red star cluster eastward of the southern nucleus 
(WS95-80) justifies it being the brightest object in the merger 
at 15/im, and how in general the 15/im emission relates to the 
expected bolometric flux. What we found is that WS95-80 is in- 
deed expected to be a bolometrically bright object, amongst the 
brightest sources in the Antennae, but not as bright as expected 
from the 15/im image. While we would expect it to contribute 
3.6% to the total flux at 15um. lMirabel et al. l dl998) measured 
it to contribute roughly 15%. Measurements concerning the 
contribution of the overlap region as a whole are more in agree- 
ment, while they measure « 50%, we would expect 43% of the 
bolometric luminosity to originate in that as 3kpc x 5kpc large 
region. 




Fig. 19. Overlay of our simulated mi w i image (slightly less 
smoothed than in Fig.l20» with radio 6cm (high-pass filtered, 
so it shows p redomin antly t he point-like sources) contours 
(from lNeff & Ulvestad 1 12000: their Fig. 4c). The radio contin- 
uum emission here consists mostly of free-free emission of gas 
photo-ionized by hot stars, therefore it is considered a good 
tracer of recent star formation. Note the excellent correspon- 
dence between the two maps. 



4. Discussion 

4.1. Age Distribution of the Young Clusters 

As shown in Section l3~5l based on Figs. [9] and ^2 about 
70% of the Kj-band detected star clusters with masses >10 5 
M0 are younger than 10 Myrs. Thus, 10 Myrs is approxi- 
mately an e-folding time for ages of the most massive clus- 
ters. Unfortunately, our information about lower mass clusters 
is severely biased by the form and sensitivity of mass-to-light 
ratio in the K-band with age. 

There are two plausible hypotheses for explaining this 
rather narrow range of ages of the Antennae clusters: Either that 
at least 70% of the clusters have been produced within the last 
10 Myrs, or that the e-folding survival time for the ensemble of 
clusters is about 10 Myrs and that either the clusters are born 
unbound or some process efficiently destroys what are other- 
wise bound clusters at birth. We note again that we can only 
firmly make this statement for clusters above »T0 5 Mq, which 
lie above the Ks-completeness limit out to ages of around 25 
Myrs. But even considering incompleteness, the number den- 
sity of older clusters below our completeness limit still decline 
and thus it is probably safe to extend the age limit for determin- 
ing the number of clusters out to an age of * 100 Myrs. 
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Fig. 20. O verlay of our sim ulated mi,„i image with CO con- 
tours (from Wilson et al. .1200 111 . The mb„i image is slightly less 
smoothed here than in Fig Rather good agreement is seen 
between the two maps in the overlap region, whereas the north- 
western region shows a significant offset between the clusters 
and the CO peaks. Note that the three regions marked in Fig. 
II 8lall lie in regions where a star cluster lies just on the edge of 
a relatively sharp edge in the CO map. The black ellipse shows 
the beam size of the CO observations. 

4.1.1. Cluster Age Distribution as Star-Formation 
History 

The changing gravitational potential in a strong merger like the 
Ante nnae can induce strong variations in th e star-formation rate 
(e.g., iMihos. Bothun. & Richstone I Il993l) . According to the 
most common models, closest passage of the two galaxies com- 
prising the Antennae happened roughly 200 Myrs ago (see e.g., 
iBarnes 119881) . At first passage pericenter, galaxy merger mod- 
els predict strong star-formation which declines and depending 
on the exact geometry and gas distribut ion, weaker peaks in the 
star- formation can o ccur (e.g., IMihos. Bothun. & Richstonel 
1 19931 Barnes 1 120041) . This enhanced star-formation tends to 
persist for a factor of 0.5 to 1 dynamical times (a "dynam- 
ical time" being the orbital time-scale of the first encounter, 
thus around 200-300 Myr) after the first passage pericenter. 
However, none of the models appear to show intense star-bursts 
that return to the star-formation rates before the strong interac- 
tion occurred. The contrast between peaks in the star-formation 
and the low points during the merger range about a factor of a 
few and always well above the rate before the strong interac- 
tion occurred. In addition, none of the peaks modelled appear 
to last less than 5% of dynamical time. 

The e-folding timescale of 10 Myrs which we estimate for 
the clusters currently observed in the Antennae is short com- 
pare d to the d ynamical time of the An tennae. In the model 
of IMihos. Bothun. & Richstone I dl993l) . the periapse of the 
Antennae occurred about 200 Myrs ago and the galaxies will 



merge in about 100 Myrs. Thus the e-folding time of the cluster 
ages represents something like less than 5% of the dynamical 
time of the merger. This is short. The disks of NGC 4039/39 
should have hardly felt a change in the potential. Nevertheless, 
clusters with identical ages are distributed throughout the face 
of the system, which is hard to explain as a "coordinated" burst, 
because the communication time between the regions is larger 
than the age spread. 

It is possible that the interaction has produced other sources 
of perturbations like spiral density waves in the disks or bars. 
In principal these could also act to induce star-formation. But 
again, the dynamical time scale for spiral density waves is sim- 
ilar to the dynamical time of the disks themselves which are 
of the order of 100 Myrs or longer. Bars can have shorter time 
scales (10s of Myrs), but would lead to star-formation which is 
largely circum-nuclear. One of the striking aspects of the star- 
formation in the Antennae is that it occurs throughout the disks. 
Moreover, we do not observe a large range of ages across the 
face of the Antennae as might be expected by gravitational per- 
turbations. Since the timescales are long, we would expect the 
age range to be a significant fraction of these timescales. While 
we do observe that the clusters are younger in the "overlap" 
region, this is only at the level of a few Myrs. Given the dispar- 
ities in the various dynamical timescales, the relatively short 
e-folding time of the cluster ages, and the narrow range of ages 
across the Antennae, it appears unlikely that the narrow range 
of ages reflects the star-formation history within the Antennae. 

4.1 .2. Cluster Age Distribution as Cluster Destruction 
Time-scale 

On the other hand, there are mechanisms which could lead to 
the destruction of clusters on a variety of time scales. The pos- 
sible mechanisms include (roughly in order of increasing time 
scale): clusters being born un- or marginally bound, rapid loss 
of natal material, mass loss through stellar evolution, stellar 
ejection during 2-body relaxation, clusters being born with a 
significant fraction of their binding energy in "hard" binaries, 
tidal interaction with nearby clusters or gravitational perturba- 
tions such as bars or spiral arms, and gravitational shocks. Of 
course, rapid destruction as constrained by the clusters in the 
Antennae may require several of these to act in unison. 

The simplest hypothesis to explain the short timescale dur- 
ing which the clusters must be disrupted would be that they are 
not gravitationally bound when they form. If true, this would 
present a very interesting quandary given that they are very 
dense stellar systems. A model for star-formation has been pro- 
posed where the driving force is supersonic turbulent conver- 
ge nce flows in the in terstellar medium (see the excellent review 
bv lMac Low & Klessen 2004; and references therein). A situ- 
ation like that, as opposed to gravitational collapse with frag- 
mentation, may lead to compact, vigorous star-formation as ob- 
served in the young compact clusters in the Antennae, but with 
the stars that formed not bound as a single entity. 

In lMengel et al. Id2002l) we estimated that the crossing time 
for massive clusters in the Antennae of tcr = R\)ilcr « 
few x 10 5 yrs. The crossing for less massive clusters than stud- 
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ied in Men gel et all J2002h will be roughly simi lar since they 
have roughly similar mean densities. Following iFall I (2004), 
if a cluster is approximately freely expanding as a results of 
catastrophic mass loss or being born unbound, the characteris- 
tic radius increases as R\iz(i) ~ R\ii(Q)( t I t cr) an d its charac- 
teristic surface density decreases as E(f) * £(0)(t/tcr)~ 2 . Thus 
after 10-100 cross times, the characteristic radius has increased 
by that amount, and the characteristic surface density would be 
a factor of 100 to 10000 lower (5-20 magnitudes). The size of 
these effects means that clusters would drop below our detec- 
tion limit in only a few crossing times. 

The magnitude of this effect mean s that it can be tested ob- 
servationally. In Menge l et al. I (|2002) we argued from an anal- 
ysis of the compactness, velocity dispersion, and IMF of some 
of the most massive young (<10 Myrs) clusters in the Antennae 
that the clusters must be relatively long lived (survive for more 
than several 100 Myrs). This conclusion was reached under the 
assumption that the clusters are initially bound and then subse- 
quent stellar evolutionary and dynamical processes are the only 
effects that lead to the dissolution of the clusters. This study 
was also limited to a small number of clusters and certainly 
could have been biased to those that are long lived (the ages 
of the clusters studied were all close to 10 Myrs). Our results 
presented here allow us to address this question with greater 
statistical power. If the cluster were unbound or dissolving due 
to internal mass loss (stellar evolutionary and dynamical), we 
would expect the cluster to expand and its concentration to go 
down with age. As argued previously, if the clusters were ini- 
tially unbound and freely expanding the concentration would 
decrease rapidly with time. In an attempt to determine whether 
the cluster size or concentration changes with age, which might 
be an indication of cluster disruption if older clusters are gen- 
erally more expanded or shallower than younger clusters, we 
measured cluster sizes and King concentration parameters for 
clusters in two of our age bins for the K-detected clusters: 
younger than 4 Myrs (7 members measured) and between 8 and 
1 1 Myrs (11 members measured). In principle, we would have 
liked more and narrower age bins and checked for an evolution 
of size with age, but there are not enough I-band bright, iso- 
lated clusters to fill these bins in statistically relevant numbers, 
therefore we used these two clearly separate age bins. 

Our result was that the spread in cluster sizes and concen- 
trations is large for both groups (in the results, we quote the 
standard deviation), but that the average sizes are smaller for 
the older group, and they tend to be more concentrated. Young 
clusters: R c // = 16 pc (+15pc), King concentrations 15-40, 
older clusters: R e ff = 6.5 pc (+5.3pc), King concentrations 
7-27. If the clusters were simply freely expanding, the effect 
should be strong and obvious. Thus free expansion seems un- 
likely. These results do not lead to a clear interpretation for 
other processes that might lead to rapid dissolution and whose 
effects are less obvious. Either, cluster evolution generally con- 
centrates the clusters, or the weakly concentrated clusters get 
disrupted, which decreases the size average for the survivors. 

For clusters which are not very concentrated after for- 
mation and/or had a shallow IMF, relaxation processes can 
lead to destruction on timescales of a few tens of Myrs (see , 
e.g., Takahashi & Porte gies Zwart 11200 0. Meng el et al. 1 12002). 



Relaxation processes could lead to rapid cluster destruction 
if the cluster is only marginally bound, as might be the case 
say for a strong convergenc e flow like formation mechanism 
iMac Low & Klessen I (y004). Given there is no clear depen- 
dence with age on any of the parameters measured, we have 
no unambiguous answer what process is dominant. Processes 
that lead to rapid dissolution short of free expansion are stel- 
lar mass loss (winds and supernova) that reduce both the stel- 
lar mass and remove the natal material (in conjunction with 
the radiation pressure), and relaxation processes (e.g., 2-body 
relaxation, interactions with hard binaries) that preferentially 
eject the low mass stars. Dissolution of clusters is enhanced for 
clusters with shallow IMFs and low stellar concentrations. So 
a plausible explanation for the generally rapid dissolution of a 
significant fraction of the clusters is that some are marginally 
bound and thus the processes such as loss of the natal material, 
stellar evolutionary-driven mass loss, and dynamical mass loss 
all become more effective and more efficient in unbinding some 
of the clusters. 

If rapid dissolution of young comp act clusters 3 — which 
is also a favoured scenario for Fall (see lFallll2004 which also 
focuses on the cluste rs in the Antennae) and for Bastian et al. 
(Bas tian et al. Il2005t which analyzes the cluster population in 
M51) — is common to all merging galaxies, then estimates 
of the fraction of star-formation in mergers (and perhaps other 
types of star-bursts) that occurs in su ch clusters is likely t o be 
severely underestimated. For example, Meur er et al. 1 lll995l) es- 
timated that about 20% of the UV light from starbursts comes 
from young compact star clusters. This would suggest that ei- 
ther there is background star-formation that does not reach the 
high stellar densities of the young compact star-clusters or that 
the background light is composed of stars from dissolved clus- 
ters. Our results indic ate that this is likely the case. Indeed, 
Tremonti et al. 1 ll200lll in a UV spectroscopic study of cluster 
and background populations in the dwarf galaxy NGC 5253 
found a relative paucity of O-stars in the background popula- 
tion compared to the compact clusters. They interpreted this 
difference as the background light (non-cluster light) being due 
to clusters that dissolved on the order of 10 Myrs. Combined 
with our direct statistical evidence fo r cluster destruction (see 
also IZhang& Falll 1 19991 lFallll2004l) provides powerful evi- 
dence for both a short lifetime of clusters and that cluster for- 
mation must be a very significant mode of star-formation in all 
star-burst galaxies. 

4.2. Nature of brightest 15/jm peak 

With its observed age of 4 Myrs and a photometric mass of 
just above 7x10 6 Mq, cluster WS95-80 is expected to have a 
bolometric observed magnitude of m,i, /=11.4 mag. This is far 
less than what is expected from the 15^/m image , where it con- 
tribut es 15% of the total flux from the merger jMirabel et al. I 
19981). This d i screp ancy is in support of a suggestion made by 
Wilson et al. (2001) for the origin of the strong 15/im emis- 

3 The first referenc e we are aware of of a dissolution timescale of 
a 10 Myr was made bv lLada & Ladal <2003l) for Galactic open clusters 
and therefore not strictly relevant for our compact, massive targets 
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sion: One of their two options was that the giant molecular 
clouds, three of which seem to be overlapping or colliding in 
this location, efficiently transport very small dust grains into 
the vicinity of the hot stars which are present in the relatively 
young (4 Myrs) star cluster. There, they are heated up to « 100 
K, which leads to strong continuum emission around 15jum 
(the spectrum is observed to be rather steeply rising in that 
region). For their other option, an even younger star cluster 
(younger than 1 Myr), we see no evidence, even in Ks-band 
or, from visual i nspec tion, in the 8yum Spitzer image presented 
bv lWang etaDj2004l) . 

We identify two other regions where we might expect a 
strongly rising continuum in the 15/im range: Roughly in the 
middle of the overlap region, and to the south-west of the north- 
ern nucleus (marked with arrows in Fig. II 81 . Both of these 
locations, like cluster WS95-80, have a 15/im flux which is 
substantially brighter than the one expected from the artificial 
bolometric image. We did not analyze the ISO data for these re- 
gions, but it could be worth while checking whether they show 
the same spectral signature like WS95-80. 

It may be significant in that context that all three regions 
which we marked in Fig. ^] are star clusters which lie just on 
the edge of a rather steep increase in CO flux. 

Concerning the coincidence and offset between CO inten- 
sity and simulated mj„; (see caption of Fig. l20t . we see this in 
the context of the age distribution in the Antennae: The overlap 
region, which hosts the youngest clusters, shows a clear cor- 
relation between CO flux and the location of the clusters (as 
traced by the bolometric luminosity), which is not the case for 
the northwestern region. There, most of the clusters are rather 
around 10 Myrs old, an age where they have blown free of their 
natal dust clouds, which can explain the observed offset. 



4.3. Extinction evolution 

In Sec. l3.6l we gave two pieces of evidence that clusters are usu- 
ally formed in high extinction regions, and have cleared away 
this natal dust cocoon by the age of 8-9 Myrs. First evidence 
is the correlation between high-extinction older clusters with 
high extinction young star clusters. This, we interpret to mean 
that a cluster of around 10 Myrs or older which suffers high ex- 
tinction is often not obscured by its own dust cocoon any more, 
but just happens to be located in a region of more recent star 
formation where the dust content is still generally high. This 
effect is not seen for clusters below 8 Myrs, and sets in around 
9 Myrs, therefore the clusters on average seem to clear away 
their natal cloud around that age. 

The same result is obtained from Fig. [H] where - despite 
large variations in Ay within each age bin - the evolution of 
extinction with age clearly declines. Also here 8-9 Myrs is ap- 
proximately the limiting age above which many clusters have 
a sufficiently low extinction that one can be sure that they have 
emerged from their natal dust cocoon. 



4.4. Impact of extinction 

As shown in the previous section, owing to extinction, a careful 
analysis of the optical data is required if sensible parameters, 
like for example total masses and thereby star formation rates, 
are to be derived from optical wavelengths alone in galaxies 
like the Antennae. 

Here, a combination of considerably deeper optical images 
and a moderate (Ay « 1 .3) average extinction means that only a 
small fraction of clusters which are detected in Ks is undetected 
in the other bands. This holds even for U band (see Table l3.6l . 

However, converting the observed cluster magnitudes into a 
quantity like mass without information on the extinction leads 
to unacceptable uncertainties. This is even the case if age infor- 
mation of the single clusters is taken into account. In that case, 
the mass, assuming zero extinction, will always be underesti- 
mated. In our case, only a quarter of the mass are derived from 
U and B observations. 

Even worse results are obtained if average ages and extinc- 
tions are assumed, then for more or less reasonable assump- 
tions for these averages, the derived total mass can easily vary 
by orders of magnitude. This behaviour will be less extreme 
in galaxies where the major star formation event is not equally 
young, because after a few hundred megayears, the luminosity 
evolution of a star cluster is much slower. 

On the bright side, we showed in the previous section that 
it is possible to obtain reasonable estimates of both, ages and 
extinction, of these simple stellar populations by fitting the 
observed broadband colours with these two free parameters. 
However, it is crucial there to have a broad wavelength base- 
line, therefore optimally one needs a detection in both, U and 
K band, or at least in one of them. 

4.5. Luminosity- and mass function 

As shown in Sec. 13.71 we constructed luminosity- and mass 
functions for our Ks-detected clusters. Even though both func- 
tions can in principle be fit by broken power laws with turnover 
cluster masses of 2.4 and 3 x 1O 5 M0, respectively, we caution 
that both function do not allow for this straightforward inter- 
pretation. 

While the luminosity function is easily constructed and also 
completeness correction is relatively easy, it is not very mean- 
ingful to compare this luminosity function to, for example, 
that of globular cluster systems, because the strong variation 
of mass-to-light ratios between the clusters building up this lu- 
minosity function means that a given luminosity cannot reli- 
ably be converted into a cluster mass. The observed "turnover 
mass", where the two power laws which were fit to the func- 
tion cross, is converted into a characteristic mass assuming the 
mass-to-light ratio of the median cluster mass in our sample, 
which is 9.2 Myrs, but particularly the high luminosity clusters 
are often rather around 6 Myrs old, where the mass-to-light ra- 
tion differs from that at 9.2 Myrs substantially. 

The mass function has a similar problem: While it is sim- 
ple to construct a mass function from those clusters which had 
ages and extinctions assigned and therefore an individual mass 
estimate, it is difficult to apply a completeness correction. Even 
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though above » 10 5 Mq, the mass function should be complete 
for all ages, the completeness limit for the median age lies 
around 3x10 4 Mq, but in the mass range between these two val- 
ues, the mass function will be affected by those clusters which 
have ages where the completeness limit lies in that mass range. 
And from the mass function strictly above 1O 5 M0, one cannot 
really reliably determine a turnover, expected to lie between 1 
and a few times 10 5 Mq- Data which do not accurately probe 
the turn-over provide very little constraint on the overall mass 
function. 

So even though we constructed the luminosity- and mass 
functions, mai nly for the sake of comparison with what o ther 
authors found dWhitmore et al. Il999tlzhang & Fall 1 19991 for 
the luminosity and mass function, respectively), and find good 
agreement at least for the luminosity function, we caution not to 
take these results too seriously. To improve the results, K-band 
observations would be required which go « 2.5 mag deeper, 
pushing the completeness limit down to around 10 4 Mq for all 
ages. And we note that without the near-infrared data, accu- 
rately accounting for the extinction is challenging, making it 
difficult to estimate robust total stellar masses for individual 
clusters. 

Even though this condition is fulfilled for t he optical HST 
data a lone (from which the mass function in IZhang & Fall I 
ill999il . which could be fit by a single power law, was con- 
structed), we think that the broader wavelength coverage ob- 
tained by including the K-band substantially improves the age- 
and extinction-, and thereby the photometric mass determina- 
tions. 

5. Summary 

To summarize, our main conclusion from the analysis of 
VLT/NTT Antennae near-infrared images, in combination with 
the HST optical images, is that there is evidence that the ma- 
jority of the bright, young cluster population will not evolve 
into globular cluster, but rather disrupt on timescales as short 
as tens of Myrs. If this was a generic property of starbursts, 
the amount of star formation in similar events would have been 
underestimated so far. 

Large uncertainties are also introduced into estimates of 
star formation rates if individual cluster ages and extinction are 
not taken into account. In detail: 

- Most Ks-band detected star clusters are young, the median 
age is w 9 Myrs, with a fairly narrow age range. Most of 
the youngest clusters (those below « 5 Myrs) are found in 
the so-called overlap region between the two galaxies. All 
other ages, also the few older clusters (up to 200 Myrs) 
are relatively evenly distributed over the two galaxies. The 
numbers of young clusters peak at two ages, 6 and 10 Myrs, 
the latter being mostly a selection effect (the clusters were 
detected in K-band). 

- The most likely explanation for the narrow range of ages is 
that clusters are destroyed or dissolve on scales of a few 10 
Myrs. An analysis of the cluster concentrations as a func- 
tion of age does not reveal a trend. Thus it is unlikely that 
the clusters are born freely expanding. More likely is that 



they are marginally bound and stellar evolutionary and dy- 
namical mass loss lead to clusters becoming quickly un- 
bound. Thus whether a cluster dissolves is likely to be de- 
pendent on its mass distribution (concentration, etc) and 
initial mass function. 

- The extinction is variable, between Ay »0 and more than 
10 mag foreground screen extinction. The average value is 
1.3 mag. It drops with increasing cluster age, which is con- 
sistent with the picture that a star cluster forms in a very 
dense and dust enshrouded environment and blows free of 
this cocoon after a few Myrs. Most clusters have reached 
this stage between the age of 8 and 9 Myrs. 

- Comparing photometry in the different bands, we note that 
information about age and extinction for individual clus- 
ters is necessary if meaningful conclusions about proper- 
ties like total stellar mass in clusters or (clustered) star for- 
mation rate are to be obtained. Not taking into account this 
information, or assumption of average values can lead to re- 
sults differing by orders of magnitudes. However, both age 
and extinction information can be obtained from broadband 
colours alone if the wavelength coverage is broad enough. 

- A lower limit to the star formation rate over the last 25 
Myrs is «20 Mq yr _1 , a relatively typical value compared 
to those for other nearby starbursts. But this is taking into 
account only the detected clusters and not corrected for any 
potential dissolution of clusters - including the latter, the 
lower limit for the SFR increases to ^50 Mq yr _1 . 

- Even though the Ks-band luminosity- and mass functions 
are best fit by broken power laws, and the turnover masses 
at face value lie around (up to factor of 2) the characteristic 
mass of globular cluster mass functions, this result should 
not be over-interpreted, due to limitations implied by sys- 
tematic uncertainties. 
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